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Abstract. Many theories of modified gravity, including the well studied Horndeski models, 
are characterized by a screening mechanism that ensures that standard gravity is recov¬ 
ered near astrophysical bodies. In a recently introduced class of gravitational theories that 
goes beyond Horndeski, it has been found that new derivative interactions lead to a partial 
breaking of the Vainshtein screening mechanism inside any gravitational source, although 
not outside. We study the impact of this new type of deviation from standard gravity on 
the density profile of a spherically symmetric matter distribution, in the nonrelativistic limit. 
For simplicity, we consider a polytropic equation of state and derive the modifications to the 
standard Lane-Emden equations. We also show the existence of a universal upper bound on 
the amplitude of this type of modified gravity, independently of the details of the equation 
of state. 
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1 Introduction 

It is one of the biggest challenges of modern cosmology to understand the physical origin of the 
cosmic acceleration of the Universe discovered in Refs. [1, 2]. It might eventually require the 
presence of a new type of matter, usually called dark energy, or alternatively a modification of 
general relativity on cosmological scales. In addition to the standard tensor modes, models 
of the latter type often involve a new scalar degree of freedom that accounts for cosmic 
acceleration. Among them, the simplest modification is described by a scalar-tensor theory 
where a single scalar degree of freedom is added (see e.g. [3]). A very general class of scalar- 
tensor actions was found in [4] under the assumption that it leads to second-order equations 
of motion for both the metric and the scalar field. These theories, usually called Horndeski 
theories, include a large number of concrete models of modified gravity as specific cases. Since 
modifications of gravity are strongly constrained by precision tests of gravity on small scales, 
such as in the Solar system, any modified gravity model involving an additional propagating 
degree of freedom must include a mechanism that suppresses any potential fifth force that 
might appear on small scales. The Vainshtein mechanism [5] represents one of the screening 
mechanisms, which is known to operate in models with nonlinear derivative interactions for 
the scalar mode. These interactions can be large in the vicinity of a gravitational source, 
leading to self-screening of the fifth force. The general theory exhibiting Vainshtein screening 
mechanism for the Horndeski Lagrangian has been studied in the literature [6-8] (see also 
[9, 10]). 

Until recently, Horndeski theories were believed to be the most general healthy scalar- 
tensor theories because of the generic occurence of Ostrogradski instabilities in systems with 
higher order equations of motion. However, it was proposed in Refs. [11, 12] that Horndeski 
theories could be further generalized without introducing dangerous instabilities, after a 
detailed study of the true propagating degrees of freedom (see also [13] for a confirmation of 
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these results). As found in Ref. [14], a specific feature of this new class of theories is that it 
exhibits a deviation from general relativity inside a region filled with matter, while the usual 
Vainshtein screening is recovered just outside. 

Considering a spherical object with a radial density profile denoted p(r), the radial 
gravitational force in the nonrelativistic regime can be expressed as 

d$ (M d 2 VW\ 
dr “ Gn V r 2 C dr 2 J 

where A4(r) = 47r f/ r' 2 p(r')dr' is the mass enclosed in the sphere of radius r. The last term 
on the right hand side represents the deviation with respect to Newton’s law. It is propor¬ 
tional to the radial derivative of the density, i.e. dp/dr and its amplitude is characterized 
by the parameter e, which can be written explicitly in terms of the functions that appear 
in the Lagrangian beyond Horndeski (See Ref. [14] and Appendix A for details). Although 
our main motivation for studying the gravitational law (1.1) comes from the models beyond 
Horndeski introduced in Refs. [11, 12], the main discussion in this paper is independent of the 
underlying gravitational theory and thus applies to any theory leading to a phenomenological 
modification of this type. 

The expected consequence of a deviation of Newton’s law is a modification of the inter¬ 
nal structure of any astrophysical body. This implies that modified gravity can in principle 
be constrained from observations of astrophysical objects whose internal physics is well un¬ 
derstood, such as stars. The goal of the present work is to study, without entering into 
the astrophysical details, how the density profile of a spherical object would be qualitatively 
modified by a gravitational force of the form (1.1). 

For simplicity, we will assume the equation of state to be polytropic. In standard 
gravity, the density profile for polytropic equations of state is determined by solutions of the 
so-called Lane-Emden equation [15]. With the new gravitational interaction given in (1.1), 
we obtain a modified Lane-Emden equation, which can be solved numerically. We also show 
that, independently of the equation of state for matter, there is a universal upper bound on 
the value of e in order to obtain physically sensible solutions. 

This paper is organized as follows: In section 2, we derive the modified Lane-Emden 
equation. We then show the existence of a critical value for e. In the subsequent section, we 
perform the numerical integration of the modified Lane-Emden equation. We also construct 
an exact analytic solution for a specific value of the polytropic index. Section 4 is devoted 
to our conclusions. 
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2 Density profile in modified gravity 


2.1 Modified Lane-Emden equation 

In this section, we consider a static, spherically symmetric distribution of matter as a simple 
toy model for an astrophysical object such as a non-relativistic star. Although this model 
is too simple to be directly confronted with observational data, it enables us to capture the 
essential modifications of the stellar structure and to provide simple estimates of the novel 
effects due to the breaking of the Vainshtein screening. For a static, spherically symmetric 
non-relativistic source, the hydrostatic equation reads 


dP d$ (M d 2 M\ 

i P i P I 9 ^ 2 ~ ) 5 

dr dr \ r z d r z J 


( 2 . 1 ) 
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where P is the pressure and we have used the modified gravitational equation (1.1). Ex¬ 
cept the latter, all other equations governing the matter distribution are supposed to be 
unchanged. Note that in our treatment, e is a priori assumed to take any value, both posi¬ 
tive or negative, (see Appendix A for the possible values of e in the context of the theories 
beyond Horndeski.) 

In order to get rid of the integral function A4(r), one can multiply (2.1) by r 2 /p and 
then take its derivative, which yields 

12 

= 0. (2.2) 


_d 

dr 


r 2 dP\ 2 

- 7 — + 47rG N r 

p dr J 


,, R v R d P 2 d P 

(1 — be) p — b e r- - er 

dr dr^ 


To obtain a closed equation and be able to solve it explicitly for the radial profile, one needs 
to specify the equation of state for matter. Here, we simply assume a polytropic equation of 
state of the form 


P = Kp 1 


(2.3) 


where K and n are positive constants. 

Following the standard procedure, we introduce the following dimensionless variable £ 
and function % such that 


£ = -, p=p c [xm n , 


where p c denotes the energy density at the center and 


r c = 


(n + 1 )Kp, 


-i+4 


47tGn 


(2.4) 


(2.5) 


Using these quantities and combining eqs. (2.1)—(2.5), after straightforward calculations, we 
obtain a modified Lane-Emden equation, given by 


^_d 

£2d£ 




= -x 


( 2 . 6 ) 


For e = 0, one recognizes the standard Lane-Emden equation. 

The density profile can be obtained by solving this equation with the following boundary 
conditions for x(£) at the center of a star £ = 0 : y(0) = 1, (d\/d£) |^ =f) = 0. Once 
the solution of this equation is computed, the radius and mass of the object, respectively 
denoted R and M, are determined by the first zero £i of the function y, i.e. x(£i) = 0; 
through R = r c £\ and the integration of the density up to the radius, 


M = 4-7T r^Pc 



£ 2 [x(£Fd£. 


(2.7) 


Note that R/r c and M/Anr^p c are uniquely determined for a given value of e. Therefore, the 

3 —n 

usual scaling M oc R 1 ~ n is unchanged, although the proportionality coefficient is modified, 
as it depends on e. Using eq. (2.6) and the definition of £i, the integration for the mass (2.7) 
can be analytically performed and gives 


M 


A 3 c 2 d X 
—4vrr c p c £, — 

A 3 c 2 d X 
-47rr c p c £, — 


f=Ci 



n > 1 

- £ f?) , 


n = 1 


( 2 . 8 ) 
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For n > 1, this result is independent of e and coincides with the expression obtained from the 
standard Lane-Emden equation. By contrast, the mass diverges for n < 1, which indicates 
that the modification of gravity leads to a dramatic change in this case. In the following, we 
will assume n > 1. 


2.2 Universal bound on e 

As we will see explicitly in the next section by solving numerically the modified Lane-Emden 
equation, it turns out that there exists a critical value for the parameter e beyond which one 
cannot find a meaningful profile. As we now show, this critical value for e can be derived 
from very general arguments that are independent of the details of the equation of state. 

Any nonsingular solution is characterized by density and pressure profiles with finite 
values and vanishing first derivatives at the center. This means that one can expand p(r ) 
and P(r) near the center r = 0, as 


1 r 2 

P = Pc + - 2 P2-^ + 


1 r 2 

P = Pc+ 2 P ' 2 lP + 


(2.9) 


Inserting these expansions into (2.2), one obtains, from the term at lowest order in r, the 
relation 


P2 = 


4irG N p 2 R 2 
-3-(1 - 6e) 


( 2 . 10 ) 


For any physically reasonable equation of state, one expects the density and pressure to 
decrease when going away from the center, which means in particular P 2 < 0. The above 
relation shows that this is possible only if e < 1/6. 

Physically, this can be understood by noting that gravity becomes effectively repulsive 
at the center of the object when 1 — 6e < 0. Indeed, near the center, we have 



d 2 M M 
-~ 6 —— . 

dr 2 r 2 


( 2 . 11 ) 


Substituting these into the right hand side of (1.1), one finds that the gravitational force is 
proportional to (1 — 6e) and thus changes sign when e crosses the critical value 1/6. 

Given a physically reasonable equation of state for which the pressure increases with 
the density, the above result that the pressure increases with r when e > 1/6 implies that 
the density also increases with r near the center. In fact, one can also show that the density 
gradient cannot vanish at any higher radius (see Appendix B for details). Consequently, this 
puzzling behavior is not confined to the core and continues for higher radii. We thus conclude 
that one cannot construct a physically sensible profile when e > 1/6. 


3 Solutions for the density profile 
3.1 Numerical solutions 

Let us first discuss the asymptotic behavior of x around £ = 0 to avoid the numerical 
instability at £ = 0. Expanding \ around the center of the core with the appropriate boundary 
condition, namely y(£) = )T) m =o c m £ m with co = l,ci = 0, we find that the following 
coefficients solve the equation (2.6), 

Tl 

C 2 = -x + e, c 3 = 0 , c 4 = — (1 - 6e) (1-20e) , c 5 = 0 , • • • . (3.1) 

6 12U 
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We note that the nonvanishing coefficients of c m approach to zero as e —>• 1/6. Hence we 
should take care of the convergence when performing the numerical integration. 

We now present the results obtained by solving the modified Lane-Emden equation (2.6) 
numerically. We consider the polytropes with the index n = 3 and n = 1, which give models 
for main sequence stars and neutron stars, respectively. For these values of ra, we investigate 
how the new gravitational force changes the density profile of a star for a given value of e. 




Figure 1. Numerical results for the profile %(£) with n = 3 (Left) and the one with n = 1 (Right). 
Different color of the lines corresponds to the different value of e as is shown in the left panncl. 


In Fig. 1, we plot the solutions of the modified Lane-Emden equation with the polytropic 
index n = 3 (Left) and n = 1 (Right) for several values of e. We can find two features from 
these figures: 

• For e larger than the critical value 1/6, the density never approaches to zero. This 
implies that a physically sensible profile cannot be obtained, as expected from the 
discussion in the previous section where it was found that the force is repulsive in this 
case. 

• By contrast, for e below the critical value 1/6, the force is always attractive (see Fig. 
2). One can see that the profile near the core tends to be steeper for smaller e, while it 
becomes flatter when moving further out. This indicates that gravity becomes stronger 
in the inner region but weaker in the outer region. 




Figure 2. The gravitational force and its sign (Left) and relative difference in the magnitude of the 
force (Right) for n = 3 and various values of e, with the same colors as in Fig. 1. In particular, e = 0 
is the light blue line. Note that because of eqs. (2.1), (2.3) and (2.4), the magnitude of the force is 
directly related to the steepness of the profiles in Fig. 1. 


- 5 - 












Fig. 2 confirms that when e is negative, the magnitude of the force 


|-fgrav| — Gj<fP 

is larger than in standard gravity near the center, as expected since the boundary conditions 
imply d 2 A4/dr 2 > 0. Further out, however, d 2 .A4/d?’ 2 becomes negative as it also contains 
dp/dr which is negative. Consequently, the magnitude of the force is reduced with respect to 
the standard situation, thus leading to a larger radius of the object. In the case 0 < e < 1/6, 
one finds the opposite effects. 

In Fig. 3, we plot = R/r c and M / (drtr/.p c ) as a function of e, for three equations of 
state (n = 1.5, 2 and 3. We observe that the dimensionless radius R/r c always decreases as e 
increases. By contrast, the variation of the dimensionless mass M /(47rr/? p c ) depends on the 
equation of state. 


M 

r 2 


d Z M 


— €- 


dr 2 


(3.2) 




Figure 3. Numerical results for the dimensionless radius = R/r c (Left) and the dimensionless 
mass M/(47rr^/9 c ) (Right) as functions of e. e = 0 corresponds to the case with standard gravity. As 
illustrative values of the polytropic index, we choose n = 1.5 (green), n = 2 (blue), and n = 3 (red). 


3.2 Analytic Solution 

In standard gravity, it is well known that the Lane-Emden equation can be solved analytically 
in the case n = 1. The corresponding analytic solution is given by 


x(0 


sin£ 

T" 


(3.3) 


In this section, we show that analytic solutions can be found also for the modified Lane- 
Emden equation (2.6) in the case n = 1. The analytic solutions are helpful to point out 
the existence of the critical point more explicitly as well as to check the numerical solutions 
obtained before. 

For n = 1, by changing the variable from £ to x defined by x = e£ 2 , the equation (2.6) 
can be rewritten as 


x(l 


d 2 y /3 7 \ dx 

x ) d ?- + (2 - 2 / s 


25 4 + e\ 

16 ~ 16e ) 
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X = 0. 


(3.4) 














We see that eq. (3.4) takes the form of Euler’s hypergeometric differential equation and the 
solution satisfying the boundary conditions is given by 


x(£) = 2^1 


1 /4 + e 5 1 /4 + e 

4V~’ 4 4 V e 



( 3 - 5 ) 


where 2-^1 (a, b, c; x) is the hyper geometric function. 

As is expected, these solutions satisfy eq. (3.1) around the center. Furthermore, we 
can see explicitly that 2-^1 (0,5/2, 3/2, £ 2 /6) = 1 holds identically, which confirms that the 
value e = 1/6 corresponds to a threshold beyond which one cannot find a physically sensible 
profile. We have also checked that the profile x(0 obtained numerically for n = 1 coincides 
with the analytical expression (3.5). 


4 Summary 

In this paper, we have investigated some consequences of a partial breaking of the Vainshtein 
screening mechanism inside an astrophysical object, which has been shown to arise in theories 
beyond Horndeski. Assuming a polytropic equation of state, we have derived a diffential 
equation for the density, which generalizes the Lane-Emden equation obtained in standard 
gravity [eq. (2.6)]. By solving this equation numerically, we have obtained the modified 
density profile, and, consequently, the modified mass and radius of a spherical object. The 
modifications depend on the amplitude of the new gravitational effect, characterized by a 
parameter e. We have found the existence of a critical value for e, below which viable profiles 
do not exist. We have argued that this is also true for any physically reasonable equation of 
state, i.e. the equation of state where the density increases with pressure. This result can be 
physically understood from the fact that gravity effectively becomes repulsive for e > 1/6. 
This is why a stable distribution of matter cannot be reached. 

In our analysis, we have considered a very simple set-up: a static, spherically symmetric 
and non-relativistic object with a polytropic equation of state. To make detailed comparisons 
with observations, and thus obtain precise constraints on the allowed amplitude of this new 
effect, it would be necessary to refine our analysis by including a realistic equation of state, 
as well as rotation. Relativistic corrections should also be taken into account for compact 
stars. 


Note added: While this paper was in preparation, Ref. [18] appeared, which also investigates 
the consequences of a partial breaking of the Vainshtein screening mechanism inside an 
astrophysical object. We note that they restricted their attention to a limited region of 
parameter space corresponding to e < 0 in our notation. 
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A Expression of the e parameter 
A.l Connection to the EFT parameters 

Let us be more explicit about the form of e in eq. (1.1). For concreteness, we study the 
quartic Lagrangian 

L 4 = G^,X)R + G iX {^X) [(D ^) 2 - - — t’ X) , 

(A.l) 

with the shorthand notation 

= V p,( f> , fypLV — V V V ^ (f>. 

The gradient of the gravitational potential is given in this case by [14] 

dd> (M_ a* d 2 M \ 

dr N \ r 2 4 ttA4 2 jG n S dr 2 J ’ 

where Mp\ is a mass scale for the rnetic perturbations introduced in [14] 
and eq. (1.1) gives the expression of e in terms of the model parameters: 


4ttM^G n E ' 

The various terms are those defined in [14] and take the specific form of 

^a*ElF 4) (8nG N )- 1 =2G 4 -8X(G 4 x+XG 4 xx)-^X 2 (5F 4 + 2XF 4X ) (A.5) 

A 15 

E = Q (4ai«2 — 20 : 1 a* + Qu) — 2 Ta\ , Mp^Q = 2 (G 4 — 2AG 4 jy) , Mp-^J- = 2G 4 (A.6) 
and 

= Gix + ZXGixx + A(5F 4 + 2XF iX ), 

= G^x + AF 4 , 

= Gixx + 2F 4 + XF iX ■ 

Mpi to be, 

Mpj = M 2 = 2G 4 - 4 AG 4 x — 4A 2 F 4 , (A.10) 

which corresponds to the mass scale M in [12]. It canonically normalizes the tensor pertur¬ 
bations and should be always positive to avoid ghost instability. 

It is also convenient to reexpress e in the language of the so-called EFT formalism, 
which has been applied to Horndeski theories in [16] and to their extensions in [12] (see also 


(A.7) 

(A.8) 

(A.9) 


Mpi 

A 3 ai 

Mpi 

AV“ 2 

V 

A6 

Here, we set the energy scale 


(A.2) 

(A.3) 

. Comparing this 

(A.4) 
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[17] for a recent review). In terms of the parameters used in [12], defined as 


as — 


4X 

MS 


[4 XF 4 + 2 X 2 F ax + G ax + 2XG 4X x] , 


4X 

«T = + XF4] > 


a H = 


4X 2 

JlS 


■F 4 , 


we find that the functions introduced earlier take the form 


J 7 = 1 + OiT , Q = 4 + OLH ) 


_ Off 

M*A 3 ) a * ~ 4 ’ 


Gn = [8nM*(l + as)} , 


and 


X \ an — cub 

Oil = --- 


M* A 3 
X 




OiT 


M*A 3 ) “ 2 4 


V — 


OLH — &T — OLB 


M*A 3 1 

Finally, one can express eq. (A.4) as 


8 


a 


e = 


H 


OLB — OiT — OIb( 1 + OIT ) 


(A.ll) 
(A.12) 
(A. 13) 

(A-14) 
(A.15) 

(A.16) 

(A.17) 

(A.18) 


A.2 Viable parameter region 

In the previous subsection, we expressed the e parameter in terms of four EFT parameters. 
To avoid several inconsistencies, these parameters cannot be arbitrary. First, because the 
gravitational force around a star is attractive, the effective gravitational constant G x should 
be positive, 


(87tGn) 1 = M* (1 + ag) > 0, (A.19) 

or 

/* = (SttGnM* 2 )- 1 = 1 + a B > 0 . (A.20) 

Here, the dimensionless parameter /* can be understood as the strength of the graviton 
interaction with respect to that of the gravitational force around a star. In addition, we 
should require that the speed of sound of the tensor perturbation should be positive, 

Ct = 1 + OiT > 0 , (A.21) 


to avoid the gradient instability. In terms of these positive quantities, the e parameter is 
expressed as, 


e 


a 

1 + OiH 


2 

H 



(A.22) 
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Therefore, e can be either positive or negative depending on the value of a#. 1 
When G A = Mpj/2 and F .4 = const, the e parameter is simply expressed as, 


e 


X 2 F a 


(A.23) 


In the G 3 -galileon theory considered in [18], where F A = — 1/A 6 , it is always negative. The e 
parameter corresponds to 4T in [18]. 2 


B No physically sensible profile for e > 1/6 


For a physically reasonable equation of state where the density increases with pressure, 
eq. ( 2 . 10 ) implies that the density is also a monotonically increasing function of r near the 
center when e > 1/6. Here, we show that the density continues to increase away from the 
center and then no physically sensible profile is obtained. Let us assume there exists a radius 
r* where the density turns to decrease, meaning (dp/dr) (r*) = 0. Since p is assumed to have 
been strictly increasing for r < r*, the following inequality holds 

47T q 

M{r*) K —p(r*)rl, (B.l) 

because p(r) has the maximum at the outermost radius r*. Moreover, by using the relation 
between mass and density, we also obtain, 


d 2 M_ 

dr 2 


(r*) = 4-7T 


2 p(r*)r, 


+ ^(r*)r 2 


8 vrp(r*)r* . 


(B.2) 


Inserting this, along with eq. (B.l), in the hydrostatic equation (2.1) yields 

> -L^A (1 -6 £)> 0. (B.3) 

which also implies (dp/dr) (r*) > 0 under our assumption on the equation of state. This 
contradicts the assumption (dp/dr)(r*) = 0. Therefore, unless the relation between pressure 
and density is an exotic one, there is no possibility of having an astrophysical body when 
e > 1 / 6 . 


dP 

dr 
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